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y—* ' Abstract 

o . 

CN ' In this paper, we obtained a new class of axially symmetric cos- 

mological mesonic stiff fluid models in the context of Lyra's geome- 

I try. Expressions for the energy, pressure and the massless scalar held 

^ ' are derived by considering the time dependent displacement field. We 

■ found that the mesonic scalar field depends on only t coordinate. Some 

^ I physical properties of the obtained models are discussed. 
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^ : 1 Introduction 

> : 

2 ■ After Einstein proposed his theory of general relativity he succeeded in ge- 

if^ , ometrizing the phenomena of gravitation by expressing gravitational in terms 

^ J I of metric tensor gij. This idea of geometrizing gravitation inspired many 

G^ I physicists to generalize the theory in order to incorporate electrodynamics as 

^ ■ a purely geometrical construct. One of the first attempts in this direction was 

O ■ made in 1918 by Weyl [1] who suggested a theory based on a generalization of 

>■ , Riemannian geometry, by formulating a new kind of gauge theory involving 

^ i metric tensor, to geometrize gravitation and electromagnetism. This the- 

5^ I ory was criticized due to non-integrability of length of vector under parallel 

displacement 0. Also as pointed out by Einstein that this theory implies 
that frequency of spectral lines emitted by atoms would not remain con- 
stant but would depend on their past histories, which is in contradiction to 
observed uniformity of their properties [2]. Nevertheless, completely apart 
from these criticisms, Weyl's geometry provides an interesting example of 
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^In the theory of general relativity if a vector undergoes a parallel displacement, its 
direction may change, but not its length. While in Weyl's geometry not only the direction 
but the length may change and depends on the path between two points, which show that 
the length is not integrable. 
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non-Riemannian connections. Folland |3] gave a global formulation of Weyl 
manifolds clarifying considerably many of Weyl's basic ideas thereby. 

In 1951 Lyra [1] proposed a modification of Riemannian geometry by in- 
troducing a gauge function into the structure less manifold which is in close 
resemblance to Weyl's geometry [5]. This modification was to overcome the 
problems appeared in Weyl's geometry and is more in keeping with the spirit 
of Einstein's principle of geometrization, since both the scalar and tensor 
fields have more or less intrinsic geometrical significance. In this way Rie- 
mannian geometry was given a new modification and the modified geometry 
was named as Lyra's geometry. For physical motivation of Lyra's geometry 
we refer to the literature [6]- [8]. However, in contrast to Weyl's geometry, 
in Lyra's geometry the connection is metric preserving and length transfers 
are integrable as in Riemannian geometry. Subsequently, Sen [9] and Sen 
and Dunn [2] proposed a new scalar tensor theory of gravitation. They con- 
structed an analog of Einstein field equation based in Lyra's geometry, see 
equation (12. 6p . Sen [9] found that static model with finite density in Lyra's 
geometry is similar to the static Einstein model, but a signification differ- 
ences was that the model exhibited red shift. Halford [10] pointed out that 
the constant displacement vector field in Lyra's geometry plays the role 
of a cosmological constant in the normal general relativistic treatment. Hal- 
ford [11] showed that the scalar tensor treatment based in Lyra's geometry 
predicts the same effects, within observational limits, as in Einstein theory. 
Several attempts have been made to cast the scalar tensor theory of gravita- 
tion in wider geometrical context [T2] . 

Many Authors [13] have studied cosmological models based on Lyra's ge- 
ometry with a constant displacement field vector in the time-direction. Singh 
and his collaborators [H] have studied Bianchi type I, HI, Kantowski-Sachs 
and new class of models with a time dependent displacement field. They 
have made a comparative study of Robert son- Walker models with a constant 
deceleration parameter in Einstein's theory with a cosmological terms and in 
the cosmological theory based on Lyra's geometry. 

Recently, several authors [I5]- [21] studied cosmological models based on 
Lyra's geometry in various contexts. With these motivations, in this paper, 
we obtained exact solutions of the field equations for mesonic stiff fluid models 
in axially symmetric space-times within the frame work of Lyra's geometry 
for time varying displacement field vector. 

Axially symmetric cosmological models have been studied in both Rie- 
mannian and Lyra geometries. In context of general relativity theory, by 
adopting the comoving coordinate system, these models with string dust 
cloud source are studied by Bhattacharaya and Karade [22]. They shown 
that some of these models are singular free even at an initial epoch. In the 
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context of Lyra's geometry these models are studied in the presence of cosmic 
string source and thick domain walls [23] and in the presence of perfect fluid 
distribution [21]. 

This paper is organized as follows: The metric and field equations are pre- 
sented in section 2. Section 3 deals with solving the field equations. Finally, 
in section 4, concluding remarks are given. 



2 Fundamental Concepts and Field Equations 

Consider the axially symmetric metric [22] in the form 

ds' = df - A\t){dx' + f\x)d<P') - B\t)dz\ (2.1) 

with the convention x° = t,x^ = x, x^ = (^, x'^ = 2; and A and B are func- 
tions of t only while / is a function of the coordinate x only. 

The volume element of the model (12.11) is given by 

V = = A^fB (2.2) 

The four-acceleration vector, the rotation, the expansion scalar and the shear 
scalar characterizing the four velocity vector field, m", respectively, have the 
usual definitions as given by Raychaudhuri [25] 

e = w^, 

where ^ 

o-jj = U(^i-j) + U{iUj) - -&{gij + UiUj). 

In view of the metric (12. ip . the four-acceleration vector, the rotation, the ex- 
pansion scalar and the shear scalar given by (I2.3p can be written in a co moving 
coordinates system as 



Ui = 0, 

= Kii(i)'+5(lr+^). 



2A , B (2.4) 
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The non vanishing components of the shear tensorcTjj are 



an = AilQA-A), 
CT22 = AfilQA-A), 

an = -|e. 

The field equations in normal gauge for Lyra's geometry as obtained by 
Sen [9] (in gravitational units c = 8ttG = 1) read as 

1 3 3 

^ij - 2^9ij = -Tij - + -gij(l)a(p", (2.6) 

the left hand side is the usual Einstein tensor, whereas (pi is a time-like 
displacement field vector defined by 

0, = (0,0,0, A(t)), 

and Tij is the energy momentum tensor corresponding to perfect fiuid and 
massless mesonic scalar field and is given by 

Tij = (p + p)u,u, - pg,j + ViV, - ^gijVkV'''. (2.7) 

Here p is the pressure, p the energy density and Ui the four velocity vector 
satisfying the relation in co-moving coordinate system 

giju'u^ = l, u' = Ui = (1,0,0,0). 

However, V is the massless scalar field and we assume it to be a function of t 
and X coordinates. The Scalar field V is governed by the Klein-Gordan wave 
equation 

g''v,j = 0. 

For the line element (12.11) . the field equations (12.61) with equation (12. 7p 
lead to the following system of equations 

A B AB 3^^ , 1-2, , , 

+ -A2 = -(p + -y2) + --, (2.9) 



A B AB 4 1 ' 1A^' 
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A,^ 2AB f" 3 , 1 • 

2A B 

P+(P + P)(^ + ;g)=0, (2.12) 

V - —V" + {— + —)V - = (2.13) 

Here the over heat dot denotes differentiation with respect to t and over head 
prime denotes differentiation with respect to x- 
From equations (12. 8p and (12. 9p . we get 

V = 0. (2.14) 

Consequently, the mesonic scalar field does not exit in the direction of x- 

The functional dependence of the metric together with equations (12. 9p 
and (I2.10p . using equation f l2.14p . imply that 

f" 

— = k^, k'^ = constant. (2.15) 

If /c = 0, then the solution of this differential equation is f{x) = kiX + ^2, 
ki and k2 are constants of integration. Without loss of generality, we choose 
ki = 1 and k2 = 0- Thus we shall have 

fix) = X- (2.16) 

In the case /(x) = x the field equations (I2.8p - (l2.13p . using equation 
(12141) . reduce to 

A B AB 3^^ , 1 • , , 

^ + (^f+jA' = -(p+5n (2^18) 

2A B 

P+iP + P)i^ + ^)=0, (2.20) 
9 A B 

i>+(— + -)y = 0. (2.21) 
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We can easily find from (12.211) that 



V = ^, (2.22) 

wliere 0) is a constant of integration. 
From equations (I2.17P and fl2.18p . we get 

B AB A A,^ 

We assume A to be some arbitrary function of B, say 

A = ^{B). (2.24) 

So equation (I2.23P becomes 



dip 

~^ H^A — 

(i-1) 



where ip^ = j^. Equation (I2.25P results in the following possibilities. 



(i-2) 



(i-3) 



i? = 0, and ^ + (^^-1^ = 0, (2.27) 

%p ip B%p 



B = 0. (2.28) 



3 Solutions of Field Equations 

We have only five highly non-linear field equations fl2.17p - fl2.21l) in sex un- 
knowns , A, B,p, p, V and A. In order to obtain its exact solution, we assume 
one more physically reasonable condition amongst these variables. We con- 
sider here the effective "stiff fluid" distribution, that is, a perfect fluid with 
the equation of state: 

p = p. (3.1) 



6 



The equation of state (13.11) was apparently first proposed by Zeldovicli [26] . 
It sliould liave applied in the early Universe, the justification being the ob- 
servation that with (13. ip the velocity of sound equals the velocity of light, so 
no material in this Universe could be more stiff. 

Using the condition (13. ip in equation (I2.20p and by integrating, we get 

P = P=JI^^ (3-2) 

where m{y^ 0) is a constant of integration. 
Case (i-1): 

From equation (I2.26p . using the first equation in the second equation, then 
by integrating, we get 

tP = ciB + C2, (3.3) 

where Ci{^ 0) and C2 are integration constants. Using this result in the first 
equation in (12.261) . we get C2 = 0. So that equation (13. 3p becomes 

ij = CiB. (3.4) 

Using this result in equation (I2.25p . we have 

A = ciB. (3.5) 

Now using this equation and the condition (13.11) . in equations (l2.17p -f l2TT9|) . 

we get 

IB , 5 , o 3 , o 3 , , 



These equations yield 



5 = (at + b)3, (3. 



where 0) and h are constants of integration. 

According to equations (13. 5 p and (13.81) the line element (12. ip can be written 
in the following form 

ds^ = dt^ - {at + h) § {cldx^ + clx^d(i? + dz^), (3.9) 

Physical properties of the model 
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Using equations (13 .Sp and fl3.8p in equations fl2.19p - fl2.2ip . take into ac- 
count (13. ip . the expressions for density p, pressure p, massless scalar field V 
and displacement field A are given by 



m 

(at + 6)2' ^ 4 



P=P = TZr—^^-: = TI' 



which shows that p and p are not singular, 



n 



V = — ^ log(at + 6) + C3, 

2 _ C4 _ 2(2m + n2) 

~ (at + 6)2' ''"""g 3^1 ■ 

It is observed, from equations (13. 5p and (13. Sp that A{t) and B{t) can be sin- 
gular only for t — )■ oo. Thus the line element (13. 9p is singular free even at 
t = 0. 

For the line element (13.91) . using equations (12. 2p . (12. 4p and (12. 5p . we have 
the following physical properties: 
The volume element is 

V = clxiat + b). 

This equation shows that the volume increases as the time increases, that is, 
the model (13. 9p is expanding with time. 

The expansion scalar, which determines the volume behavior of the fluid, is 
given by 

at + 

The only non- vanishing component of the shear tensor, aij, is 

2a 



a 44 



3{at + b)' 
Hence, the shear scalar a is given by 

^3(at + 6)^ ■ 

Since \imt^oo{^) 7^ 0, then the model (13. 9p does not approach isotropy 
for large value of t. Also the model does not admit acceleration and rotation, 
since Ui = and Uij = . 

Case (i-2): 
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The first equation in f l2.27p gives 

B = bit + b2, (3.10) 

where ^1(7^ 0) and 62 are constants of integration. Using this result in the 
second equation in (12.271) and take into account (I2.24p . we get 

= bsihf + 2b2t) + bi, 

where 63(7^ 0) and 64 are constants of integration. This equation can be 
written in the form 

A = ib,B^ + b,)K (3.11) 

where the constant 65 depends on the constants bi and f)2- 

Using equations (13.101) and (13.111) in equations (I2.17p - (I2.21I) . take into account 

the condition p = p, we get 

^ = ^=(MT6^[^^^^^^ + ^^)' + '^]'' 



where b^ is the constant of integration, 



A 



2 2{2blbl{bit + b2f - 2m - v?) Sbjb-s 



3{b,{b,t + b2Y + b^Yibit + 62)2 3(63(61^ + 62)2 + 65) ■ 

In this case the hne element (12. ip takes the following form 

ds^ = df - (&3(&it + ^2)' + h){dx^ + x^dcf") - {b,t + b2ydz\ (3.12) 

We shall now give the expression for kinematic quantities. A straightfor- 
ward calculation leads to the expressions for element volume V, expansion 
scalar 6 and shear tensor aij of model (I3.12p are given, respectively, by 

V = x{bit + b2){bs{b,t + b2f + b,), 
which shows that the model is expanding with time, 

2bMbit + b2) , b. 



e 



0-11 



ihib,t + b2y + b,) (M + &2)' 

biibsibit + b2y + b^) bMbit + b2) 



3(61^ + 62) 
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0^22 = X^O-ii, 



2bMbit + b2f 2bi{ht + b2 

(^33 ^ 



£744 



b3{bit + b2y + b, 3 
46163(61^ + 62) 26i 



3(63(61^ + 62)2 + 65) 3(61^ + 62)' 

and other components of the shear tensor, aij, being zero. Hence 



a' 



1/ 6f6i(6it + 62)2 . 61 .2 26^63 
9V (61^ + 62)2 + 65)2^ ^6it + 62^ 63(61^ + 62)2 + 65 



Moreover, this model represents non-rotating and has vanishing acceleration. 
Since limt_5.oo(^) 7^ 0, then the model fl3.12p do not approach isotropy for 
large value of t. 
Case (i-3): 

Equation f l2.28p can be easily integrated to give 

B = e, (3.13) 

where 0) is the constant of integration. Using this equation in equations 
fl2.17l) - fl2.19l) and take into account the condition fl3.2p .we have 

A={iit + i2)K (3.14) 

where £i{^ 0) and £2 are constants of integration. Using equations fl3.13p 
and (13.141) in equations fl2.17p -f l2^2T|) . take into account the condition p = p, 
we get 

m 

p = p 



P{ilt + i2 



-2) 

n 



\2- 



y = — iog(4t + ^2) + 4, 

where £3 is the constant of integration, 

^ 2(2m + n2) 



(£1 + ^2)2' ^ V ^ 3^2 • 

In this case the line element (12. ip can be written in the form 

ds'^ = dt^ - {£it + £2){dx^ + x^d(l)'^) - £^dz\ (3.15) 

From equation (I3.14p . one can observed that A{t) is singular only when 
t — )■ 00. Consequently, the line element (I3.15P is singular free even t = 0. 
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For the line element (13.151) . the volume element V and the kinematics prop- 
erties (acceleration -Uj, rotation Uij, expansion scalar 0, shear tensor aij and 
shear scalar a) respectively found to have the following expressions: 



which shows that the model is expanding with time, 

w = 0, 



Uij = 0, 



e 



o 

(722 = xVii, 



0"33 



0-44 



36(4t + £2) ■ 

As in the above two cases, limt^oo(§) 7^ 0, the model (I3.15P do not approach 
isotropy for large value of t. 



4 Conclusions 

This paper deals with axially symmetric space-time in the presence of mesonic 
stiff fluid distribution within the framework of Lyra's geometry for time de- 
pendent displacement field. We have presented a new class of exact solutions 
of Einstein's field equations for this space-time. The obtained models rep- 
resent shearing, non-rotating and expanding with time t. Moreover, these 
models are singular free even at the initial epoch t = and have vanishing 
accelerations. For all models, we found also that \im.t-^^{^) ^ 0, this means 
that they are not approach isotropy for large time t. We found also that 
the mesonic scalar field in axially symmetric space-time exists only in the 
t-direction. 



11 



References 

[1] Weyl, H., Sitzungsber, preuss. Akad.Wlss. 465 (1918). 

[2] Sen, D. K. and Dunn, K. A., J. Math. Phys. 12, 578 (1971). 

[3] Folland, G., J. Diff. Geom. 4, 145 (1970). 

[4] Lyra G., Math. Z. 54, 52 (1951). 

[5] Scheibe, E., Math. Z. 57, 65 (1952). 

[6] Sen, D. K., Can. Math. Bull. 3. 255 (1960). 

[7] Manoukian, E. B., Phys. Rev. D5, 2915 (1972). 

[8] Hudgin, R. H., J. Math. Phys. 14, 1794 (1973). 

[9] Sen, D. K., Z. Phys. 149, 311 (1957). 
[10] Halford, W. D., Austr. J. Phys. 23, 863 (1970). 
[11] Halford, W. D., J. Math. Phys. 13, 1699 (1972). 

[12] Brans, C. and Dicke, R. H., Phys. Rev. 124, 925 (1961); Wheeler, J. A., 
"Geometrodynamics" (Academic press. Inc. New York) (1962); Peters, 
P., J. Math. Phys. 10, 1029 (1969); Ross, D. K., Phys. Rev. D5, 284 
(1972); Dunn, K. A., J. Math. Phys. 15, 2229 (1974); Dunn, K. A., 
Tensor (N. S.) 29, 214 (1975). 

[13] Sen, D. K. and Vanstone, J. R., J. Math. Phys. 13, 990 (1972); Bhamra, 
K. S., Aust. J. Phys. 27, 541 (1974); Karade, T. M. and Borikar, S. M., 
Gen. Relativ. Gravit. 9, 431 (1978); Kalyanshetti, S. B. and Waghmode, 
B. B., Gen. Relativ. Gravit. 14, 823 (1982); Besham, A., Astrophys. 
Space Sci. 127, 189 (1986); Reddy, D. R. K. and Venkateswarlu, R. 
Astrophys. Space Sci. 136, 191 (1987); Soleng, H. H., Gen. Relativ. 
Gravit. 19, 1213 (1987); 

[14] Singh, T. and Singh, G. P., J. Math. Phys. 32, 2456 (1991) ;I1. Nuovo Ci- 
mento B106, 617 (1991); Int. J. Theor. Phys. 31, 1433 (1992); Fortschr. 
Phys. 41, 737 (199); Singh, G. P. and Desikan, K., Pramana-J Phys., 
49, 205 (1997). 

[15] Paradhan, A., J. Math. Phys. 50, 022501 (2009) 



12 



[16] Casana, R., Melo, C. and Pimentel, B., Astrophys. Space Sci. 305, 125 
(2006). 

[17] Rahaman, F, Bhui, B. and Bag, G., Astrophys. Space Sci. 295, 507 

(2005) . 

[18] Bali, R. and Chandani, N. K. J. Math. Phys. 49, 032502 (2008) 

[19] Kumar, S. and Singh C. P., Int. J. Mod. Phys., A23, 813 (2008). 

[20] Rao, V. U. M., Vinutha, T. and Santhi, M. V., Astrophys. Space Sci. 
314, 213 (2008). 

[21] Singh, J. K., Astrophys. Space Sci. 314, 361 (2008). 

[22] Bhattacharaya, S and Karade, T. M., Astrophys. Space Sci. 202, 69 
(1993). 

[23] Reddy, D. R. K. and Rao, M. V. S, Astrophys. Space Sci. 302, 157 

(2006) . 

[24] Rao, V. U. M. and Vinutha, T., Astrophys. Space Sci. 319, 161 (2009). 

[25] Raychaudhuri, A. K., "Theoretical Cosmology" (Clarendon, Oxford) 
(1979). 

[26] Zeldovich, Ya. B., Sov. Phys. JETP 14, 1143 (1962). 



13 



